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Abstract 



We investigate whether flavor changing neutrino interactions (FCNIs) can 
be sufficiently large to provide a viable solution to the atmospheric neutrino 
problem. Effective operators induced by heavy boson exchange that allow 
for flavor changing neutrino scattering off quarks or electrons are related by 
an SU{2)l rotation to operators that induce anomalous tau decays. Since 
SU{2)l violation is small for New Physics at or above the weak scale, one 
can use the upper bounds on lepton flavor violating tau decays or on lepton 
universality violation to put severe, model-independent bounds on the relevant 
non-standard neutrino interactions. Also Z-induced flavor changing neutral 
currents, due to heavy singlet neutrinos, are too small to be relevant for the 
atmospheric neutrino anomaly. We conclude that the FCNI solution to the 
atmospheric neutrino problem is ruled out. 
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I. INTRODUCTION 



Several atmospheric neutrino (AN) experiments [jT]-Q observe an anomalous ratio Uf./ fjj, 
in the atmospheric neutrino flux. This long-standing AN problem has been conflrmed by the 
recent Super-Kamiokande high-statistics observations |Q, which give strong evidence that 
the standard model (SM) description of the neutrino sector is incomplete. The standard 
solution to the AN anomaly in terms of neutrino oscillations requires that neutrinos are 
massive, and that there is mixing in the lepton sector. Then, z/^ oscillations can 

explain the atmospheric neutrino data provided that the relevant mass-squared difference 
is A^Tv ~ 10~^ eV^ and the muon and tau neutrinos have large vacuum mixing angles, 
sin26'^Ar ~ 1. 

Recently, an alternative solution, where the AN anomaly is induced by non-standard 
neutrino interactions has been proposed In this scenario the neutrinos are assumed 

to be massless, but they are subject to non-standard interactions. For neutrino propagation 
in matter, flavor changing neutrino interactions (FCNIs) induce an off-diagonal term in the 
effective neutrino mass matrix, while non-universal flavor diagonal interactions generate the 
required splitting between the diagonal terms. A priori such a scenario is well motivated, 
since many extension of the SM predict new neutrino interactions. Moreover, it is well known 
that one cannot explain the atmospheric 0, solar and LSND [Tl| neutrino anomalies 



with three light neutrinos. Thus, rather than ignoring one of the results or introducing a 
forth, light sterile- neutrino |T^, it is interesting to investigate whether FCNIs can explain 



any of the three neutrino anomalies . 

The two effective parameters that describe the non-standard interactions of and Vt 
are '^'^'^ ^ 



14,15 



el = ^ and e'l = "^"^ ^ , (1.1) 

where G^^^,^ (a, j3 = fj,,T and f = u,d, e) denotes the effective coupling of the four fermion 
operator 

Ot^{y:;up)Uf). (1.2) 

The Lorentz structure of Ol depends on the New Physics that induces this operator. Op- 
erators which involve only left-handed neutrinos (and which conserve total lepton number) 
can be decomposed into a {V — A) ® {V — A) and a. (V — A) ® (V + A) component. (Any 
single New Physics contribution that is induced by chiral interactions yields only one of 
these two components.) It is, however, important to note that only the vector part of the 
background fermion current affects the neutrino propagation for an unpolarized medium at 
rest [0. Hence only the {V — A) ® {V) part of Ol is relevant for neutrino oscillations in 
normal matter. One mechanism to induce such operators is due to the exchange of heavy 
bosons that appear in various extensions of the standard model. An alternative mechanism 
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arises when extending the fermionic sector of the standard model and is due to Z-induced 
flavor-changing neutral currents (FCNCs) [|l^,|l^]. 



Recent analyses of non-standard neutrino interactions as a possible solution to the 
atmospheric neutrino data suggest that FCNIs can provide a good flt to the data provided 
that 

el,e'l> 0.1 or el,e'l>0.3. (1.3) 

In Ref. only new interactions involving the d quark were considered. Since the earth 
is electrically neutral and its neutron to proton ratio is close to unity, we conclude that the 
required values for el and e'l are similar for q = d,u, while those for and e'l are larger by 
a factor of three. 



The authors of Ref. |1T9| , |20| argue that such a scenario does not lead to a good description 
of the data. In this paper we do not try to resolve this issue, but investigate whether the 
lower bounds on el and e'l in ( p.. 3D are at all phenomenologically viable. The authors of 
Ref. have estimated the upper bound < 0.1 — 1.0 from the low energy z/^ neutral 
current cross-section, concluding that FCNIs could be relevant for the AN anomaly. In 



Ref. 1^ speciflc models that could give el and e'l as large as in (|1.3|) were discussed. We 



argue, however, that model-independently the upper bounds from related, charged lepton 
decay data imply that el or e'l can be at most at the one-percent level. Thus we conclude 
that FCNIs do not play a signiflcant role for the atmospheric neutrino problem. 

In Section ^ we investigate in a model-independent framework the constraints on FCNIs 
that are induced by heavy boson exchange. In most cases the upper bounds on lepton flavor 
violating tau decays, in particular t' ^ fi^ M (M = 7i^,p^,ri) and —>■ fi~ e~ , imply 
stringent constraints on el that are inconsistent with (|1.3| ). In the remaining cases severe 
constraints on e'l are derived using bounds on lepton universality violation. In Section ^ 
we show that also Z-induced FCNCs, that arise due to heavy singlet neutrinos, cannot be 
large enough to explain the AN anomaly. We conclude in Section 



II. FLAVOR CHANGING NEUTRINO INTERACTIONS INDUCED BY HEAVY 

BOSON EXCHANGE 

A. Formalism 

The analysis of FCNIs that could be relevant for the AN problem is similar to the 
discussion in Ref. [0, where the possibility that FCNIs explain the LSND results [0] was 
ruled out. In general, the presence of a heavy boson B that couples to fermion bilinears Bij 
with the trilinear couplings Xij, where i,j = 1, 2, 3 refer to fermion generations, gives rise to 
the four fermion operator Bj-Bki with the effective coupling 
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at energies well below the boson mass Mg. Thus, in terms of the trilinear coupling Aq,/ 
that describes the coupling of some heavy boson B to z/^ (a = /i, r) and a charged fermion 
f = u,d,e the effective parameters in (LI) are given by 



The crucial point of our analysis is the following: Since the SM neutrinos are components 
of SU{2)l doublets, the same trilinear couplings X^f that give rise to non-zero el or e'l also 
induce other four-fermion operators. These operators involve the SU{2)l partners of the 
neutrinos, i.e. the charged leptons, and can be used to constrain the relevant couplings. 
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Tab. 1: Lepton-Lepton Bilinears 
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Tab. 2: Lepton-Quark Bilinears 
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Tab. 3: Quark-Quark Bilinears 
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In order to obtain a complete list of these operators we note that Lorentz invariance 
implies that any fermionic bilinear Bij can couple to either a scalar (S) or a vector (V) 
boson. If the two fermions of the bilinear have the same (opposite) chirality they require 
scalar (vector) couplings. To form a gauge invariant trilinear coupling, the boson B must 
have opposite hypercharge Y and transform in the appropriate representation of SU{2)l 
and SU{3)c- Since the SM neutrinos only appear in doublets of SU{2)l and since all right- 
handed (left-handed) charged fermions transform as SU{2)l singlets (doublets), it follows 
that any boson B that couples to the fermionic bilinear can only be a singlet (s), a doublet {d) 
or a triplet (t) of SU{2)l. All relevant bilinears containing only leptons are listed in Tab. 1, 
and those that are built from a lepton doublet and a quark are listed in Tab. 2. In Tab. 3 we 
list the relevant diquark bilinears, namely, those that can couple to LL. Here Q and L denote 
the left-handed SM quark and lepton doublets, and e^j, ur, dn refer to the right-handed SM 
singlets. Some of these couplings appear in well known extensions of the standard model. 
For example, in supersymmetric models without i?-parity (SUSY l^p) p2|, fermion bilinears 
can couple to left-handed sleptons {L'^), right-handed sleptons (-^r), left-handed squarks {Q'^) 
and right-handed down squarks ((i/j), as indicated in the third column of the tables. An 
example for a scalar triplet is the in left-right symmetric models (LRSMs) p3| . 

In general any two bilinears appearing in Tab. 1-3 that couple to the same boson can 
be combined to a four fermion interaction with effective coupling as given in Eq. ( |2.1j ). In 
order to generate a non-zero el or e'{, in Eq. (|2.2|) at least one of the bilinears has to contain 



a lepton doublet L. Clearly, four fermion operators that are the product of one bilinear and 
its hermitian conjugate can be constructed. If the two bilinears have the same (different) 
flavor structure the resulting operator will conserve (violate) lepton flavor. In addition, the 
(LL) bilinear can couple to i^eR) or to any of the quark-quark bilinears in Tab. 3 inducing 
four- Fermi interactions of the form (LL)(e^eK), {LL){QQ), {LL){urUr) and {LL){djidR). 

Note that scalar or vector fields that transform identically under the unbroken SU{3)c^ 
U (1) EM symmetry can mix. If this mixing is between a doublet and a singlet or a triplet the 
resulting operators violate total lepton number and are not relevant for our analysis. In the 
case of singlet-triplet mixing no new operators are generated. Therefore this kind of mixing 
does not affect our conclusions and we neglect it. 

To demonstrate how SU (2) l related processes can be used to constrain the parameters 
el or e'l, let us consider for example the bilinear L/r (/ = e, m, d) that couples via a scalar 
doublet to its hermitian conjugate fnL. In terms of the component fields the effective 
interaction is 

where ia = fJ^LyTi for a = ^,t. X^f is the trilinear coupling of Lafn to the scalar doublet 
and Ml 2 denote the masses of its SU{2)l components. The important point is that the 
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scalar doublet exchange not only gives rise to the four-Fermi operator Ol in ( |1.2| ) (with 
(y — A) ® {V + A) structure), but also produces the SU{2)l related operator 



0{^{Lip) iff) 



(2.4) 



which has the same Lorentz structure as Ol, but where the neutrinos are replaced by their 
charged lepton partners. Moreover, the effective coupling of 0{, that we denote by G^^, is 



related to Gl by 



cf Ml 



(2.5) 



Constructing all the relevant four fermion operators that are induced by the couplings 
between the bilinears listed in Tab. 1-3, one finds that in general is generated together 
with Ol'. Here /' can be different from / only for interactions with quarks, that is in some 
cases Ol (O^) is generated together with O^ (C^")- The leptonic operator is always 
generated together with O^ unless the interaction is mediated by an intermediate scalar 
SU (2) L singlet that couples to 



(2.6) 



where la = fJ'L, tl for a = ^,t, with the elementary coupling \ae- The coupling of (LaLe) 
to {LpL(.)\ that is mediated by a scalar singlet of mass M yields the effective interactions 



(2.8) 



M2 
2M2 



where we used a Fierz transformation and the identity A'^'y^B'^ = —B'-y^A to obtain 
One can see that in this case Of, is generated together with three more operators that have 
the same effective coupling (up to a sign). However, unlike for the case of intermediate 
doublets or triplets, all these operators involve two charged leptons and two neutrinos. 



B. Experimental constraints 



There is no experimental evidence for any non-vanishing . Therefore, whenever 0( 
is generated together with Ol, one can use the upper bounds on to derive constraints 



on Gl ^ . 

fi~ e"*" e" 



The most stringent constraint on is due to the upper bound on r 



123,121 



BR(r- 



/i e 



e-) < 1.7 X 10" 



(2.9) 
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Normalizing the above bound to the measured rate of the related lepton flavor conserving 
decay, BR(r" 



jj, U^Ur 



0.17 we obtain 

< 3.1 X 10"^ G, 



G 



(2.10) 



To constrain we may use the upper bounds on various semi-hadronic tau decays that 



violate lepton flavor |24,25 
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(2.11) 
(2.12) 
(2.13) 
(2.14) 

Let us flrst consider the tau decays into 7r° and p^. Since these mesons belong to an isospin 
triplet we can use the isospin symmetry to normalize the above bounds ( p.ll| ) and ( |2.12| ) 
by the measured rates of related lepton flavor conserving decays. Using BR(r~ VtT^^) = 

0.22 |2|J2|] we obtain 



0.11 H and BR(r 



< 8.5 X lO"'' Gf and G^^(p) < 7.5 x lO"'' Gp ■ 



(2.15) 



Since the tt (p) is a pseudoscalar (vector) meson its decay probes the axial-vector (vector) 
part of the quark current. 

In general, any semi-hadronic operator Oj can be decomposed into an / = and an / = 1 
isospin component. Only the effective coupling of the latter can be constrained by the upper 
bounds on the decays into flnal states with isovector mesons, like the vr and the p. If the 
resulting operator is dominated by the / = component, the bounds in ( p.l5|) do not hold. 
But in this case we can use the upper bound on BR(r~ p" rf) in ( p.l3| ). Since the rj is an 
isosinglet, isospin symmetry is of no use for the normalization. However, we can estimate 
the proper normalization using the relation between the rj and vr hadronic matrix elements, 
which is just the ratio of the respective decay constants, /,,//7r — 1.3 [^,^. Taking into 
account the phase space effects, we obtain from ( p.l3|) that 



GIM < 1.2 X IQ-^Gp 



(2.16) 



Since the is a pseudoscalar meson its decay probes the axial-vector part of the / = 
component of the quark current, while the neutrino propagation is only affected by the 
vector part. As we have already mentioned, for any single chiral New Physics contribution 
the vector and axial- vector parts have the same magnitude and we can use ( p^.l6[ ) to constrain 
the isosinglet component of O^. In case there are several contributions, whose axial- vector 
parts cancel each other (a vector singlet V(l, 1, 0) that couples to all the diquark singlets of 
Tab. 3 with the same strength would lead to such a scenario), the / = component could still 
be constrained by the upper bound on BR(r^ jx^ tt+tt^) in (12141 ). While the calculation 
of the rate is uncertain due to our ignorance of the spectra and the decay constants of the 
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isosinglet scalar resonances, we expect that the normahzation will be similar to that of the 
71, p and T] discussed before. Finally we note that the decay r~ ^ ji" uj would be ideal to 
constrain the / = vector part, but at present no upper bound on its rate is available. 

While one can always fine-tune some parameters in order to avoid our bounds, our basic 
assumption is that this is not the case. Thus from ( |2.10D , ( p.l5[ ) and ( ^.161 ) we conclude that 



the effective coupling could be at most at the one-percent level. 

We still have to discuss the case of the intermediate scalar singlet that couples to {LL)s 
inducing the effective interactions in ( p.8|) . Since Ol is not generated together with 
but only with operators that involve two charged leptons and two neutrinos, the decay 
fi~ e~ is of no use to constrain the effective couplings. However, since the effective 
operators in ( p.8| ) contain only left-handed fermions rather strong bounds can be derived on 
the flavor diagonal terms using lepton universality. The reason is that the corresponding 
interactions are identical to the SM ones and have to be added coherently. 



Setting a = (3 = jjL,T the last term in ( p.8|) induces additional contributions to 
^hVa'yl with the effective coupling 



I A, 



|2 



Gl ^ = . (2.17) 

These new contributions violate lepton universality and lead to a deviation of the parameter 



^ iVr(/i — > e z/^z/e) Gf 

from unity. Here N denotes a normalization factor, which is just the ratio of the above two 
rates in the SM such that Rr/fj, = 1 if G^uaUa ~ ^- ■'■^ approximation we assume that 
^ <^ Of (a = t). From the most recent experimental data |27,25] it follows that 



Rr/^, = 1.0008 ± 0.0030 , (2.19) 

implying that 

e'l = ~ ^'^'^^ < 3.8 X 10-3 . (2.20) 

Gf 

Here we used that Ol has the same effective coupling as the related operator that induces 
the new contribution to ia T^- 

Finally we remark that we can use lepton universality violation not only to constrain 
the interactions induced by an intermediate singlet, but our argument holds also whenever 
an SU{2)l related operator induces additional contributions to the SM weak interactions. 
The bounds on lepton universality violation in semi-hadronic processes are of similar 
order as the bounds for the leptonic processes that appear in the definition of -Rr/^t in ( P-18D . 



Consequently analogous arguments as those leading to the upper bound on e'^ in ( |2.20| ) can 
be used to constrain e'l- Since all involved fermions have to be left-handed this only applies 
for intermediate singlets or triplets of SU{2)l. For the triplets the effective couplings of the 
relevant operators may differ due to SU{2)l breaking effects, which we will study next. 
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C. Constraining SU{2)l breaking effects 



If SU{2)l breaking effects are negligible then G^^,^^ is equal to G{^^. Comparing the 
experimental bounds pTTq ), (gJ^, and with (|Op, we find that in the SU{2)l 



symmetric limit the new neutrino interactions that we considered cannot have a significant 
contribution to the AN anomaly. 

The excellent agreement between the SM predictions and the electroweak precision data 
implies that SU{2)l breaking effects cannot be large. To show that they cannot sufficiently 
weaken the upper bounds on el and e^-^ to be consistent with ( |1.3D , we recall from Eq. ( p.5|) 
that in general the ratio of the couplings, C^^j.^/G;^^, is given by ratio Mf/M^. Here Mi 
and M2 are the masses of the particles belonging to the SU{2)l multiplet that mediate 
the processes described by G;^^ and Gl^j^_^, respectively. If Mi 7^ M2 this multiplet will 
contribute to the oblique parameters S', U and, most importantly, T . Then we can use 



a fit to the precision data to determine the maximally allowed ratio (Mi/M2)^ax- 

We use the program GAPP by J. Erler to calculate the SM predictions. For the 
latest precision data from the Z-pole [0, the VF-boson and top quark masses [pT| , p2| , deep 
inelastic scattering [Q, neutrino-electron scattering [Q] and atomic parity violation |p5|| . 



we obtain the best fit values of the oblique parameters: 

5 = -0.07 ±0.11, T = -0.10 ±0.14, [/ = 0.11 ±0.15. (2.21) 

We calculate the contributions to S", T, and U from the various scalar representations in 
Tab. 1 and Tab. 2 and determine the best fit to the data at each value of the mass splitting. 
The best fit Higgs mass Mh varies with the splitting, and we limit the Higgs mass to the 
range 95 GeV < Mh < 1 TeV. Constructing a function we determine the upper bound on 
the mass splitting between the different members of a multiplet at a given confidence level 
(CL). The individual 90% CL bounds on (Mi/M2)^ are given in the last column of Tab. 1 
and Tab. 2. (Note that the limit on {M1/M2Y is stronger if the lightest mass is heavier. 
From the Z- width measurement, the lightest mass must be heavier than Mz/2, and the 
bounds presented in the tables are derived for the case where the lighter mass is 50 GeV.) 

We did not calculate the bounds for the vector multiplets. Since vector bosons give in 
general larger contributions, we expect the bounds in the vector cases to be as good or better 
than the corresponding bounds for the scalar multiplets. Thus, for the vector multiplets a 
rather conservative upper bound is (Mi/M2)^ < 7. 

Hence, even the maximal possible SU{2) l breaking effects could weaken the bounds we 
derived only by a factor of a few and e^, cannot exceed the few-percent level. We learn 
that also the flavor changing neutrino scattering / z^r / induced by heavy bosons, that 
are doublets or triplets of SU{2)l, cannot significantly contribute to the AN anomaly. 
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III. Z-INDUCED FLAVOR CHANGING NEUTRAL CURRENTS 



In the previous section we considered models where the only modification to the neutrino 
sector is due to new interactions mediated by heavy bosons. In this section we study the 
opposite scenario, where new fermions are added, but no extra bosons beyond the SM ones 
are needed. As an example we consider Z-induced FCNCs that arise when introducing a 
heavy sterile neutrino. Such SM gauge singlets, which appear in many extensions of the 
standard model, are frequently employed to explain the smallness of the neutrino masses 
via the see-saw mechanism |^ . 



The basic idea for Z-induced FCNCs is that if a neutrino interaction eigenstate is a linear 
combination of light and heavy mass eigenstates then the effective low-energy interaction 
eigenstates, that consist only of light mass eigenstates, are not orthogonal to each other p? 



Thus the couplings to the Z-boson (and in fact also to the M^-boson) have to be modified 
slightly, implying that also the effective Hamiltonian that describes the neutrino propagation 
in matter has to be changed. 




In Ref. WE\ a general discussion of Z-induced FCNCs and their impact on neutrino oscil- 
lations has been presented. In the context of the AN problem we illustrate the mechanism 
by considering a simple framework where besides the SM neutrinos and Ur there is only 
one new gauge singlet z/g. (For simplicity we assume that the does not play an important 
role here.) These interaction eigenstates are connected to the mass eigenstates by a unitary 
t r ansf ormation 

Url Ur2 f/.J U2 , (3.1) 
Usi Us2 UshI \l^h) 

where vi and z/2 denote the light and Vh the heavy mass eigenstates. The neutrinos that are 
produced in low-energy charged-current interactions together with charged leptons /i and r 
are 

/,,-P\ /TT TT \ /,, \ 

(3.2) 

i.e. we have projected and Vr onto the Vi — V2 plane. Since the mixing matrix appearing 
in ( p.2| ) is only a submatrix of the unitary matrix in ( p.l| ), and z/^ are not orthogonal to 
each other 

) = u;,Uri + u;,Ur2 = -u;,UrH (3.3) 

and also not properly normalized. Consequently these states do not provide a proper basis 
for the neutrino oscillation formalism. 

The description of neutrino oscillation in the presence of heavy gauge singlets has been 





worked out in Ref. [|T^]. The main result is that the effective non-unitary mixing induces a 



flavor off-diagonal contribution in the matter- induced neutrino potential VpcNC- The effect 



10 



is proportional to the neutron density and its size is characterized by the ratio between 
Vfcnc and the standard (flavor diagonal) neutral current (NC) potential Vnc- 



4 = "^^ w;„u.,.\ (3.4) 

VNC 



The approximation refers to the simple example we discussed previously. It reveals that the 
effect is in general small, since it is proportional to the components of the known neutrinos 
along the heavy mass-eigenstates, which cannot be large. The Z-induced FCNCs cannot be 
constrained by the S't/(2)L-related charged lepton decay, that we used before in the context 
of FCNIs due to heavy particle exchange, but one can obtain a stringent constraint on 
from a global fit using lepton universality, CKM unitarity, and the measured Z invisible 
decay width |^. The updated analysis in |TB[ yields the conservative bounds (at 90% CL) 



I f/^;, r < 0.0096 and |f/^;,r < 0.016 . (3.5) 

We conclude that the parameter in ( p.4| ) cannot exceed the few-percent level. Therefore 
from ( [1.3|) it follows that Z-induced FCNC effects are too small to be relevant for the AN 
problem. 



IV. CONCLUSIONS 

Extensions of the standard model in general do not conserve lepton flavor and there- 
fore provide an alternative mechanism for neutrino flavor conversion that may show up in 
neutrino oscillation experiments. While such a scenario where flavor changing neutrino in- 
teractions (FCNIs) explain one of the three neutrino anomalies is a priori well motivated. 



one has to check carefully whether these solutions are phenomenologically viable. In |T3 
it was shown that FCNIs cannot be large enough to explain the LSND anomalies. In this 
paper we argue that it is also very unlikely that the AN anomaly is due to FCNIs. Both 
analysis rely on three facts: 

• The neutrino flavor changing four fermion operator that is induced by the exchange of 
a heavy boson is related by an SU {2)l rotation to other operators that violate lepton 
flavor. 

• The strength of these related operators is severely constrained by the upper bounds 
on lepton number violating processes. 

• High precision measurements imply that the violation of the SU{2)l symmetry is not 
large for new physics at or above the weak scale. Consequently the upper bounds on 
the operators that induce FCNIs are of the same order as those of the SU{2)l related 
operators. 
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The first point follows immediately from the fact that the SM neutrinos appear in SU{2) l 
doublets. Using the upper bounds on ^ ji~ M (M = n^,p^,r]) and r~ — * fi~ e~ 
we constrain, in a model-independent way, the strength of the flavor changing neutrino 
scattering reaction f ^ f be at most at the one-percent level (compared with Gp)- 
For the unique case of an intermediate scalar singlet we derive a severe constraint on the non- 
universal flavor diagonal neutrino interactions using the upper bound on lepton universality 
violation. 

The constraints we obtained for the parameters that describe the new neutrino inter- 
actions are not consistent with the values that are required to explain the AN anomaly in 
terms of FCNIs. Thus we conclude that such a solution is ruled out. One could evade 
our bounds by fine-tuning several new physics contributions, but we do not consider such a 
scenario as very attractive. Nevertheless, we would like to stress that ultimately any alter- 
native explanation [39| of the AN anomaly should be tested by the experimental data itself. 
For the time being "standard" neutrino oscillations with massive neutrinos remain the most 
plausible and elegant solution. 
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